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Although there is not a complete “proof” of the second law of thermo¬ 
dynamics based on microscopic dynamics, two properties of Hamiltonian 
systems have been used to prove the impossibility of work extraction from a 
single thermal reservoir: Liouville’s theorem and the adiabatic invariance of 
the volume enclosed by an energy shell. In this paper we analyze these two 
properties in the Szilard engine and other systems related with the Maxwell 
demon. In particular, we recall that the enclosed volume is no longer an 
adiabatic invariant in non ergodic systems and explore the consequences of 
this on the second law. 


1 Introduction 

Since the second half of the last century it is accepted that the second law of thermo¬ 
dynamics — and consequently the Carnot principle — has a probabilistic nature and is 
only valid in average. A Brownian particle in a field, for instance, can gain potential 
energy from the surrounding fluid converting heat into work for short periods of time. 
Already in 1860 Maxwell realized that these fluctuations can be used to beat the sec¬ 
ond law and illustrated this idea with a celebrated gedanken experiment: the Maxwell 
demon . 

The Maxwell demon was the starting point for a fruitful exploration of the relation 
between entropy and information, which is rooted in the foundations of statistical me¬ 
chanics and thermodynamics. Szilard devised a stylized version of the demon, the Szilard 
engine^ in which this relation could be quantified: one bit of Information (the outcome 
of a yes/no experiment) reduces the entropy by /Mn2 and, for Instance, allows one to 
extract an energy A:Tln2 from a single thermal bath In a cyclic process. 

After Maxwell and Szilard, much effort was devoted to assess the entroplc cost of 
measurement that could restore the validity of the second law In the presence of feed¬ 
back. However, In the 1970’s, Landauer and Bennett Introduced a novel and surprising 
element: an entroplc cost In the erasure of the Information that the demon gathers In 
the measurement. Landauer showed that the erasure or, more precisely, the Isothermal 
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restoring of a memory to a reference state, dissipates heat, increasing the entropy of 
the surrounding thermal reservoir [^. Bennett went further and applied the Landauer 
principle to the Szilard engine, proving that for some setups the measurement can be 
realized at zero entropy cost and the second law is restored by the cost of erasure . 

In the last years, the thermodynamics of information has experimented a new de¬ 
velopment with the application of non-equilibrium fluctuation theorems and stochastic 
thermodynamics |4]-[^ and also due to the experimental verification of the Landauer 
principle and several experimental realizations of the Szilard engine [9-11 . 

More recently, two atypical examples of Szilard engines were introduced [12, 13 . Un¬ 
like the classical Szilard engine, which operates with a feedback loop and in contact 
with a heat bath, these systems can be cooled down during a cyclic process without 
feedback. They are called microcanonical Szilard engines, because they work if the sys¬ 
tem is prepared in a microcanonical initial state and is thermally isolated during the 
process. Despite these differences, they share a fundamental feature with the original 
Szilard engine, namely, at some stage of the cycle they undergo a symmetry breaking. 
In other words, the ergodicity of the system is lost implying that the phase space of the 
system gets split into disconnected regions. The fact that symmetry breaking is the key 
ingredient of the thermal, feedback driven Szilard engine was already pointed out in 
and explored further in [IT] . 

From the phase space point of view, two fundamental properties of Hamiltonian sys¬ 
tems forbid work extraction from a single heat source during a cyclic process: i) the 
conservation of phase space volume, known as Liouville^s theorem^ and ii) the adiabatic 
invariance of the volume enclosed by an energy shell. Therefore, during a cycle of a 
Szilard engine, at least one of these two properties must be broken. 

In this paper, we attempt an analysis of generic Maxwell demons, thermal and micro- 
canonical, from a phase space perspective. Our aim is to characterize the behavior of 
trajectory bundles in phase space under feedback and non-feedback Hamiltonian evolu¬ 
tions, in which energy shells split into disconnected regions at some stage of the process. 
We analyze the consequences of this splitting on the aforementioned properties i) and 
ii). This novel approach is particularly suited for the present special issue, since Pro¬ 
fessor Jacques Yvon contributed greatly to our understanding of the role of phase space 
dynamics in statistical mechanics through the celebrated BBGKY hierarchy. 

The paper is organized as follows. In section we review the basic properties of 
Hamiltonian dynamics in phase space with special attention to Liouville’s theorem and 
the adiabatic invariance of the volume enclosed by an energy shell. In section we 
analyze the behavior of trajectories in phase space for the original Szilard engine and 
discuss how feedback allows one to beat Liouville’s theorem and effectively contract a 
region in phase space. In section we study the consequences of the failure of the 
adiabatic invariance of the enclosed volume in non ergodic evolution and the restrictions 
that Liouville’s theorem imposes to cyclic non feedback protocols in phase space. Finally, 
in section we present our conclusions and discuss the subtle problem of the definition 
of entropy for single systems, instead of ensembles, when ergodicity is broken. 
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2 Basic phase space dynamics 


The solution of the Hamilton equations in the phase space T defines an evolution operator 

This evolution operator satisfies two restrictions which 
are characteristic of Hamiltonian dynamics: the Liouville’s theorem and the adiabatic 
invariance of the volume enclosed by an energy shell. Liouville’s theorem states that the 
volume vol(A) of a set of micro states A C T is invariant, i.e., Yol[Ut '= vol(A). 

The adiabatic invariance of the enclosed volume is more restrictive and only applies to 
ergodic systems undergoing quasistatic processes. In these processes, the system evolves 
under a time-dependent Hamiltonian iL(g,p;At), where Xt is one or several external 
parameters that change infinitely slowly in time t = [0,t], that is At ^ 0 and r ^ oc. 
For the following discussion it is useful to distinguish between an energy shell of energy 
E 

^E = {{q,p)-.H{q,p-\)=E} ( 1 ) 

and an energy layer of finite width /XE\ 

T{E-, AE) = {{q,p) : E < H{q,p- X) < E + AE). (2) 

In a quasistatic process where the system remains ergodic on every energy shell, 7 ^ is 
entirely mapped into another shell ")e^w{e)^ W{E) being the work done on the system 
along the process which generically depends on the initial energy. The work W{E) 
is deterministic because the system explores the whole energy shell at each stage of 
the process, due to quasistaticity and ergodicity together. Hence, a point {q^p) in 7 ^ 
ends the process with an energy E + W{E) lying in ")e+w{e)' Reversing the momenta 
(g,p) ^ {q^—p) and the protocol followed by the external parameter, one concludes 
that, if H{q^p] A) = H{q^ —p; A), then every point in ")eew{e) is mapped by this reverse 
process into a point in 7 ^. Therefore, we conclude that Urp["iE] — 1 eew{e)‘ Let us 
define the volume of the subset enclosed by 7 ^ as 

4>x{E) = ^ dqdpe{E - H{q,p-, A)). (3) 

Ergodicity also implies that the energy layers do not cross each other along the process, 
i.e., that E' + W{E') > E + W{E) for all E' > E. In other words, the final energy 
E + W{E) is a strictly monotonic increasing function of the initial energy E. Then 
Liouville’s theorem applied to the enclosed volume implies 

cf>x,{E) = cf>x^iE + W{E)). (4) 

In particular, for a cyclic process, Aq = A^-, VF = 0, since (l)\{E) is a one-to-one function 
of E. 

For non-equasistatic processes, the work is no longer deterministic and an energy shell 
is mapped onto a more intricate set that intersects several shells ranging from E + VFmin 
to E-hVFmax, bFrnin 2 ind VFmax being the minimum and maximum value of the work. The 
usual assumption in statistical mechanics is that the intricateness of this final set cannot 
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be described in terms of meso- or macroscopic states: any coarse graining in phase space 
effectively smears the set to cover the whole energy layer of a width Vbmax ~ bKmin- One 
can conceive other mechanisms of information loss with the same result: for instance, 
if, after completing the process, we let the system evolve for a random time with the 
time independent and ergodic Hamiltonian H{q^p] A^-), then the final set of an ensemble 
of systems starting in 7 ^ will be a layer of finite width H/max ~ bKmin- In any case, the 
effective volume of a set increases in a non-quasistatic processes. 

Summarizing, quasistatic processes map energy layers into energy layers with the same 
volume, whereas non-quasistatic processes map energy layers into intricate sets that any 
information loss or relaxation mechanism effectively convert into layers with a larger 
volume than the original one. The apparent violation of the second law by a Maxwell 
demon should be related with either Louville’s theorem or the adiabatic invariance of 
the enclosed volume 0(£’) at the level of phase space dynamics. In the next sections 
we analyze how several Maxwell demons proposed in the literature operate in the phase 
space and compromise the validity of Liouville’s theorem and the adiabatic invariance 
of the enclosed volume. 


3 The Szilard engine in phase space 

The Szilard engine consists of a single particle gas at temperature T: a particle in a 
container of volume Vq whose velocity after a collision with the walls is randomly set 
according to a Maxwell distribution at temperature T. As sketched in Fig.[^ an external 
agent or demon inserts a piston in the middle of the container, measures in which half the 
particle lies and performs a reversible expansion from a volume Vb /2 to Vq? extracting 
an amount of work 

fVo 

ITextr= / PdV = kT \n2 (5) 

JVol2 

where we have used the state equation PV = /cT, P being the average pressure exerted by 
the particle and k the Boltzmann constant. The measurement is necessary to implement 
the reversible expansion, since the external agent must exert a force identical to the 
pressure of the single-molecule gas [^. It is important to notice that the Szilard engine 
can be realized with any system that undergoes a splitting of its phase space —like 

a colloidal particle in a double well potential (TH[T^ or electrons in quantum dots —, 
using an appropriate feedback protocol. 

At the level of the phase space, the feedback control allows trajectories to cross each 
other, inducing a reduction of phase space volume. This is sketched in Fig. |^a) where 
F represents the phase space of the single molecule gas and the thermal reservoir. Any 
final microstate (qr^Pr) is reached by two different trajectories starting at U^^{qr^PT) 
and U^^{qr^Pr)^ where Ul^r is the evolution operator Urp corresponding to the protocol 
implemented when the particle lies in the left and right half of the container, respectively. 
According to Liouville’s theorem, this means that a subset in the phase space of volume 
(/)init shrinks to a final volume = 0init/2. If the final energy E+W(E) is a monotonous 
function of the initial one and the process is quasistatic, then the region enclosed by the 
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Figure 1: The Szilard engine. A demon insterts a piston in the middle of a container 
with a single-molecule gas surrounded by a thermal bath at temperature T. 
The demon measures in which half the particle lies and carries out a reversible 
expansion extracting a work /cTln2 from the thermal bath in a cycle. No¬ 
tice that the protocol followed by the demon depends on the outcome of the 
measurement. 


a) b) 

Demon state r 

R 


L 


Figure 2: Trajectories in the phase space F of the single molecule gas plus the thermal 
reservoir in the Szilard engine, a) The feedback in the cyclic protocol implies 
that two trajectories with different initial conditions trajectories cross at the 
end of the process decreasing the volume in phase space by a factor 2. b) The 
state of the demon is included in an extended phase space Ftot = F x {L, R}. 
The demon is initially in state R and the trajectories are identical as in the 
previous case a), but now the demon must be restored to state R to complete 
the cycle. 
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energy shell 7 ^, with volume </)(£’) is mapped onto the region enclosed by a shell "/e+w 
obeying 


+ W) 


<m 

2 


( 6 ) 


Taking logarithms and expanding to first order in W: 


W 

VT 


= -ln 2 


( 7 ) 


since which is the result that we already obtained applying the state 

equation of ideal gases to the single molecule gas (cfr. (© and Q, with VKextr = —W). 

Liouville’s theorem is recovered if we take into account the physical nature of the 
demon and analyze trajectories in an extended phase space Ftot = F x which 

includes the state of the demon, L or i?, depending on the outcome of the measurement. 
These are not necessarily microscopic states. They can be meso- or macroscopic states 
covering equal volumes in the phase space of the microstates of the demon (see for 
an example with a non symmetric demon). In Fig. [^b) we show trajectories in which 
the demon starts in a definite state R. The trajectories in F are identical as the ones 
discussed above, but now the final set covers, in the extended phase space Ftot, fhe same 
volume as the initial one. The total system does not complete a cycle, since the demon, 
originally in state i?, ends either in L or R. To complete the cycle one has to reset the 
demon to the original state R which implies a reduction of the phase space volume that 
has to be compensated by an increase in the phase space volume of the thermal bath. 
This increase is the heat dissipation associated to erasure, as dictated by Landauer’s 
principle. 


4 Beating the adiabatic invariance of 

In the derivation of equation @ and the discussion of the previous section, we have 
shown that the impossibility of extracting work in a cycle relies on two basic facts: 
Liouville’s theorem and the monotonicity of the final energy E + W{E) as a function 
of E. This monotonicity is warranted if the energy shells are compact connected sets 
along the whole process, i.e., if the shells consist of a single closed hypersurface in F. In 
that case, the shells cannot cross each other. On the other hand, the Szilard engine is 
an example of a quasistatic processes where energy shells split into several disconnected 
components. We will show that in those cases Liouville’s theorem is still valid but the 
enclosed volume 0 could be no longer invariant. 

Consider a processes where ergodicity is broken at some stage and energy shells 7 ^ 
split into disconnected regions i = 1 , 2 ,... (the number of regions could depend on the 
initial energy E). Even if the process is quasistatic, the work is no longer deterministic 
but can take different values Wi{E) with probability Pi{E)^ which is the probability that 
the system occupies region i in the non-ergodic stages of the process. 

We are interested on how the enclosed volume (/)Xq (E) is transformed under this non 
ergodic evolution. Let Ui{(^XQ{E)) = (^Xri^ + Wi{E)) be the transformation of the 
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Figure 3: Three examples of the transformation of the volume enclosed by energy shells. 

The initial volume (/)init = 0 Ao(^) is mapped into + Wi{E))\ 

a) corresponds to the microcanonical Szilard engine introduced by Vaikun- 
tanathan and Jarzynski with pi = pi = I (see Fig. Q; b) corresponds to 
the microcanonical Szilard engine introduced by Marathe and Parrondo [T^ 
with Pi = 1/2 (i = L^R) and pi = I (see Fig. [^; c) corresponds to a Szilard 
engine in contact with a thermal bath at temperature T with pL = pR = 1/2 
and pl — 2/3, pR — 1/3 (see Fig. [^. It is easy to check that in all cases the 
slopes verify Eq. 


enclosed volumes 0. Consider an initial layer F(E, AE) and let A(/)o = </>Xq{E + AE) — 
(pXoiE) be its volume. The volume of the set of initial microstates {q^p) in the layer that 
cross region i along the process is pi A(/)o. Now consider the time-reversed protocol with 
initial condition in the final layer r{E + Wi{E)^ AiE'). Let pi{E) the fraction of those 
states that go back to the initial layer T{E^AE). Liouville’s theorem implies: 

PiA(l)o = Pi [Ui{(l)Xo{E) + A0o) - Ui{(l)Xo{E ))]. (8) 

The limit Aipo 0 immediately yields 

which is a condition that any non-feedback evolution must fulfill. Other than that, the 
transformation of (pXoiE) can be in principle arbitrary. 

In Fig. we present three possible transformations Ui that correspond to physical 
models introduced in previous works. Figure]^ a) corresponds to the model introduced 
by Vaikuntanathan and Jarzynski in In this model, an isolated one-dimensional 

classical particle is subjected to a cyclically varying quartic potential. The protocol is 
sketched on Fig. The effect of the cycle is to interchange two energy layers (depicted 
in dark and light grey on Fig. [^. This process involves ergodicity breaking as can be 
seen on panel (d) of Fig. The work performed during this cycle is a deterministic 
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-2 -1 0 

(d) t = t/2 



(e) t = 3r/4 


(f) t = r 


Figure 4: The microcanonical Szilard engine by Vaikuntanathan and Jarzinski in 
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The model consists of a classical particle subjected to a one dimensional quartic 
potential U{q). The panels also show the evolution of regions in phase space 
formed by points with a given coordinate q (horizontal axis) and total energy 
(vertical axis). Panel (a): Initially, the potential is unimodal and the system 
is ergodic. The cycle proceeds as follows. Panels (b) and (c): a well is slowly 
created on the right part of the potential. Panel (d): a second well is slowly 
created on the left part of the potential. At that stage, the potential is bimodal 
and the ergodicity of the system is broken. Panel (e): the energy of the right 
well is slowly increased until the well disappears. Panel (f): The left well 
disappears in the same way and the potential recovers the initial shape. The 
net effect of this cycle is to swap the two regions of phase space shaded in dark 
and light grey on the figure. 


function of the initial energy. Therefore pi = 1. Since the protocol is symmetric in time, 
we also have = 1, leading to a slope equal to one on Fig. |^a), even though the curve 
is only piecewise continuous. 

Figure l^b) corresponds to a model introduced by Marathe and Parrondo [^. The 
cycle is depicted on Fig.[^ This model consists of a classical particle in a one dimensional 
box of size 2L. The cycle proceeds as follows. A potential barrier of height Vb is inserted 
in the middle of the box. Then, the right wall is moved quasistatically slowly to the 
right. Finally, the left wall of the box is moved to the position of the barrier. When 
the barrier is inserted, ergodicity is broken. If the energy of the particle is smaller than 
the height of the barrier, then the work performed depends on whether the particle is 
trapped on the left side or on the right side of the barrier. Each of these cases happen 
with equal probability pL = Pr = 1/2. In fact, any shell of energy lower than Vb gets 















Figure 5: Microcanonical Szilard engine introduced by Marat he and Parrondo in [l^ . 

The system consists of a classical particle confined in a one dimensional box 
and subjected to a potential that changes quasistatically. As in Fig. the 
panels show the potential and the evolution of regions in phase space. The 
cycle proceeds as follows. First, a potential barrier of height Vb is introduced 
in the middle of the box (panel (b)). Then, the right wall is slowly moved to 
the right (panels (b,c,d)). Finally, the left wall is slowly moved to the position 
of the barrier (panels (d,e,f)). The effect of this cycle is to transport the energy 
layers as depicted on panel (f). In particular, the layer of energy lower than 
Vb gets split into two disconnected layers, the dark grey regions in panel (f). 
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Figure 6: Process leading to Fig.j^c). A classical particle is enclosed in a closed container 
in weak contact with a heat bath. At some point, a partition is inserted 
splitting the container into two equal parts. The partition is then moved very 
slowly to the right until the volume of the left part of the container is 2/3 of 
the total volume. 


split into two disconnected shells in the dark grey area in Fig. |^f). However, in the 
reverse process the work performed is deterministic. Hence, = 1. If the initial energy 
of the particle is higher than Vg, then the work performed is deterministic in the forward 
and in the backward process, implying a slope equal to one, as seen in the corresponding 
region of Fig. [^b). 

The two previous examples are called microcanonical Szilard engines because they are 
able to extract energy from the system in a cyclic protocol, if the process starts with a 
microcanonical ensemble with energy such that < 0init- 

Our third example, corresponding to Fig. [^c), is macroscopic: a particle in contact 
with a thermal bath and following a non-feedback protocol which is a slight modification 
of the original Szilard engine. The system consists of a classical particle in a closed 
volume, in weak contact with a large Hamiltonian system (the heat bath). We perform 
the following cyclic process: We insert a rigid wall partitioning the volume accessible to 
the particle into two equal parts. Then we move this wall quasistatically to the right in 
such a way that, at the end of the process, the left part occupies 2/3 of the total volume 
and the right part 1/3 (see Fig. [^. At the end of the process, we remove the wall. 

Equation in this case reads 

Pi(l){E)=pi(l){E + Wi). (10) 


where VF^, i G {L,R} is the work performed on the system, pi is the probability that 
the particle is on the side i (left of right) of the wall, and pi is the probability that the 
particle is on the side i of wall in the reverse process. Hence, in our case, pi,—pn — 1/2, 
while pl = 2/3 and pn = 1/3. 

Equation (10) above can also be interpreted as an adiabatic invariance of the volume 
(j) corresponding to each ergodic component of the phase space. During the process, 
the system is ergodic on the part of the phase space where the particle is on the left 
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side of the moving wall. Therefore, the volume enclosed in the part of the 

energy shell where the particle is on the left of the moving wall, is invariant. Hence, 
</>imt,L(-E) = ^fin,L(-E + Wl). Moreover, (/>imt,L(-E) = PL(f>{E) and 4>Rn,L{E + TU) = 
Pi,<P{E + Wl), because for one free particle in a closed container of volume V, in weak 
contact with a heat bath, we have (piE, V) = V f{E), where / is some function. Similarly, 
0imt,R(-E) = PR(f>{E), and (f)fin,K{E + Wr) = PR<j){E + Wr). In the end: 



(t>{E + 

(t>{E + Wr) 


IHE) 

1<P{E) 



if the particle is in the left side 
if the particle is in the right side 


( 11 ) 


as depicted on Fig. [^c). 

The transformation of the phase space volume 0(£’) determines the work performed 
along the process. In the limit where the heat bath part is very large, i.e. VFl,r 
E q. 0 yields 

( 12 ) 

( 13 ) 


Wl = kTln- < 0, 


and 


WR = kTln->0, 


which is the result given by classical isothermal thermodynamics. The average work is 
(Wi) = A:Tln(3/v^), which is strictly positive. 

This argument can be extended to any process where the probabilities pi and pi do 
not depend on the energy E and there are no swaps of phase space regions. In that case, 
Eq. (© implies that the maps Ui{(l)) are linear, that is. 


0fin ^ (l)xA^ + ^ ^ 

0init 0Ao (^) Pi 


(14) 


For a large system, if Wi{E) <C E, we can expand the logarithm of the enclosed volume 
obtaining 

ln(pXr{E)-ln(pXoiE)+^^j^ = ln^ (15) 

which yields, for a cyclic process, Wi = kTlnpilpi. 

To interpret this equation for non cyclic processes, let us calculate the total enclosed 
volume assuming weak coupling between the system and the reservoir, i.e., that the 
Hamiltonian of the total system can be decomposed as X) = Hb{xb) + 

Hsys{x]X)^ where x and xb denote micro states of the system and the reservoir, re¬ 
spectively. The enclosed volume is 

(f)x{E) = j dxB J dx e [E - Hb{xb) - Hsys{x; X)] 

= J dx (Pb{E - Hsys{x; X)) 

= ME) I = ME)Zx (16) 
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where Z\ is the partition function of the system corresponding to Hamiltonian H^y^{x] A). 
The bulk entropy can be written as 


S^,^,^{E) = k\nc^x{E)^k\nME) 


F{X) 


(17) 


where E{X) = —kT\nZ\ is the free energy of the system at temperature T. Equation 


(15) for the work then reads 


Wi{E) - [F(A^) - F(Ao)] = fcTln^. 

Pi 


(18) 


This is a well known relationship in stochastic thermodynamics, first introduced for 
trajectory bundles in Notice that here we have derived this equation for single 

trajectories: Eq. (14) is valid in general for a trajectory starting with energy E and 


passing through the ergodic component i, whereas (18) can be seen as a way of rewriting 


(14) for Wi E and using the standard definitions of free energy and temperature for 


big systems. Finally, since both pi and pi are normalized, averaging (14) over pi one can 
check that {Wi) > E{Xr) — F(Ao), i.e., the only way of extracting work is by decreasing 
the free energy of the system, in conformance with the second law. 

Therefore, it is not possible to extract work in a cycle if pi and pi are constant and Ui is 
linear. We need either a swap of phase space regions, as in Fig.j^a), or energy-dependent 
probabilities Pi{E) and Pi{E)^ as in Fig. [^b). It is an open problem to reproduce this 
type of transformation in a macroscopic system, although it is not hard to find examples 
where Pi{E) depends on the energy. This occurs, for instance, when a critical point 
depends on temperature, which is the case for almost every phase transition in nature. 
Remarkably, exceptions are the symmetry breaking transitions induced by pistons or 
infinite barriers confining particles, which are the preferred ones in the Szilard engine 
literature. 


5 Conclusions: from ensembles to single systems 

Along the paper we have analyzed trajectories under feedback and non feedback qua¬ 
sistatic cyclic protocols breaking ergodicity. We have seen that it is possible to break the 
adiabatic invariance of the volume 0 (F) enclosed by an energy shell 7 ^, and to reduce 
0(F) in a quasistatic process. This can be done either by a feedback loop as described 
in section 1^ or through a clever protocol as described in section]^ In the latter case, it 
is even possible to systematically reduce the energy of a thermally isolated Hamiltonian 
system during a cyclic process, see Figs. |^a) and b). The examples known so far have 
only a few degrees of freedom, but this idea could in principle be implemented in a 
macroscopic system. Finding such a system or, alternatively, proving its impossibility, 
would help to get a better understanding of the foundations of statistical mechanics and 
of the second law. 

The loss of ergodicity poses another fundamental question regarding the definition of 
entropy. If two meso- or macroscopic phases coexist, then the behavior of a single system 
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significantly differs from the behavior of the ensemble. A single system chooses between 
one of the phases, whereas the ensemble can cover the two phases simultaneously. More¬ 
over, the ensemble is not necessarily in global equilibrium, since the populations in each 
phase depend on the whole process and not only on the value of the external param¬ 
eters at a given time [Tt] . For a single system in a microstate that lies in one of the 
disconnected energy layers, the set of available microstates is the corresponding sublayer 
and not the entire energy layer. The volume of this set can suddenly decrease due to a 
symmetry breaking, like the one resulting from the insertion of the piston in the original 
Szilard engine or the raising of the barrier in the microcanonical Szilard models described 
above. This amounts to a decrease of the effective entropy of the system which is not 
accompanied by extra work or dissipation m Id] . The decrease is counterbalanced if we 
want to drive the system back to its initial state, as it happens in the Szilard engine by 
virtue of the Landauer’s principle. However, if the cycle is not completed, the reduction 
of the phase space volume of the available set of microstates remains. We believe that 
this reduction should be taken into account when the second law is applied to single 
(isolated) systems. 

One extreme example is the whole universe. It is customary to extrapolate the sec¬ 
ond law and talk about the increase of entropy in the universe as a fundamental law. 
Nevertheless, the universe has undergone massive symmetry breaking transitions, like 
nucleogenesis. The corresponding reduction of the available set of microstates has not 
been considered up to our knowledge. Another example is the appearance of slow vari¬ 
ables like the position of a crystal, which is also the result of a symmetry breaking. 

As we have seen along the paper, the loss of ergodicity in Hamiltonian systems intro¬ 
duces novel aspects that affect the significance and scope of the second law. We have 
discussed here some of them: the evolution in phase space of thermal and microcanoni¬ 
cal Szilard engines, the possibility of extracting work from systems in a microcanonical 
state, and the problem of defining entropy for single systems instead of ensembles. All 
this indicates that symmetry breaking transitions and the loss of ergodicity play an im¬ 
portant role in the foundations of statistical mechanics that deserves further and deeper 
research. 


6 Aknowledgements 

This work has been financially supported by Grant ENFASIS (FIS2011-22644, Spanish 
Government). 


References 

[1] H.S. Leff, A.F. Rex, eds., MaxwelVs Demon: Entropy, Information, Computing 
(Princeton University Press, New Jersey, 1990) 

[2] R. Landauer, in MaxwelVs Demon: Entropy, Information, Computing (Princeton 
University Press, New Jersey, 1990) 


13 




[3] C. Bennett, Int. J. Theor. Phys. 21, 905 (1982) 

[4] T. Sagawa, M. Ueda, Information Thermodynamics: MaxwelTs Demon in Nonequi¬ 
librium Dynamics (Wiley-VCH, 2013), pp. 181-211, ISBN 9783527658701 

[5] K. Sekimoto, Stochastic Energetics^ Vol. 799 of Lect. Notes Phys. (Springer, Berlin 
Heidelberg, 2010) 

[6] B. Gaveau, L. Granger, M. Moreau, L.S. Schulman, Entropy 16, 3173 (2014) 

[7] J.M.R. Parrondo, J.M. Horowitz, T. Sagawa, Nat Phys 11, 131 (2015) 

[8] A. Berut, A. Arakelyan, A. Petrosyan, S. Ciliberto, R. Dillenschneider, E. Lutz, 
Nature 483 , 187 (2011) 

[9] S. Toyabe, T. Sagawa, M. Ueda, E. Muneyuki, M. Sano, Nature Phys. 6, 988 (2010) 

[10] J. Koski, V. Maisi, T. Sagawa, J.P. Pekola, Phys. Rev. Lett. 113, 030601 (2014) 

[11] E. Roldan, LA. Martinez, J.M.R. Parrondo, D. Petrov, Nature Phys. 10 , 457 (2014) 

[12] R. Marathe, J.M.R. Parrondo, Phys. Rev. Lett. 104 , 245704 (2010) 

[13] S. Vaikuntanathan, C. Jarzynski, Phys. Rev. E 83 , 061120 (2011) 

[14] J.M.R. Parrondo, Chaos 11 , 725 (2001) 

[15] R. Kawai, J.M.R. Parrondo, C.V. Den Broeck, Phys. Rev. Lett. 98, 80602 (2007) 

[16] T. Sagawa, M. Ueda, Phys. Rev. Lett. 102, 250602 (2009) 

[17] J.M. Horowitz, J.M.R. Parrondo, Acta. Phys. Pol. B 44, 803 (2013) 


14 



